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We perform extensive Monte Carlo simulations of the 10-state Potts model on quenched two-dimensional <1> 3 
gravity graphs to study the effect of quenched connectivity disorder on the phase transition, which is strongly first 
order on regular lattices. The numerical data provides strong evidence that, due to the quenched randomness, 
the discontinuous first-order phase transition of the pure model is softened to a continuous transition. 



1. INTRODUCTION 

Systems subject to quenched random disorder 
often show a completely different behaviour than 
the pure case. If the pure system has a continuous 
phase transition it is well known that quenched 
random disorder can drive the critical behaviour 
into a new universality class, or the transition can 
even be eliminated altogether (!]] . In the case of a 
first-order phase transition in the pure system the 
effect of quenched random disorder can also be 
very dramatic, with the possibility for a softening 
to a continuous transition. 

The paradigm for testing the latter prediction 
is the two-dimensional (2D) g-state Potts model 
which exhibits on regular lattices for q > 5 a first- 
order transition whose strength increases with q. 
In Ref. [|| the effect of quenched bond disorder 
was investigated for the 8-state model and it was 
found that the critical behaviour of the quenched 
model could be well described by the Onsager 
Ising model universality class. In Ref. [^) the ef- 
fect of quenched connectivity disorder was stud- 
ied by putting the 8-state model on 2D Poisso- 
nian random lattices showing that for this type 
of quenched disorder the transition stays first or- 
der. 

A different sort of connectivity disorder ap- 
pears in 2D gravity triangulations or their dual 
$ 3 graphs. In such models one is interested in 
the coupling of matter to 2D gravity, so the disor- 
der is annealed rather than quenched. Motivated 



by Wexler's mean field results for q = oo Potts 
models coupled to 2D gravity Q, simulations of 
the 10-state and 200-state Potts model coupled 
to 2D gravity gave convincing evidence H for a 
continuous transition. 

As the only quenched connectivity disorder seri- 
ously investigated to date, 2D Poissonian random 
lattices, showed no sign of softening for first-order 
transitions it is interesting to enquire whether the 
salient feature for the softening in the 2D gravity 
simulations is the annealed nature of the connec- 
tivity disorder or whether it is some intrinsic fea- 
tures of the graphs themselves. We discuss here 
results of a Monte Carlo (MC) study of the 10- 
state Potts model on quenched random lattices 
drawn from the equilibrium distribution of pure 
2D gravity triangulations (or, more exactly, their 
dual 4> 3 graphs). 

2. MODEL AND RESULTS 

We used the standard definition of the q-state 
Potts model, 



i, 



(1) 



where (ij) denotes the nearest- neighbour bonds 
of random $ 3 graphs with up to 10 000 sites. 
For each lattice size we generated 64 inde- 
pendent replica using the Tutte algorithm ||, 
and performed long simulations of the 10-state 
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model near the transition point using the single- 
cluster update algorithm. After thermalization 
we recorded measurements of the energy E and 
the magnetization M in 64 time-series files. The 
corresponding quantities per site are denoted in 
the following by e = E/N and m = M/N. 

Given this raw data we employed standard 
reweighting techniques to compute, e.g., the spe- 
cific heat, C«(/3) = (3 2 N({e 2 ) - (e) 2 ), for each 
replica labeled by the superindex (i), and then 
performed the replica average, C(0) = [C^ (/?)] = 
(1/64) Y.T C (4) (/3), denoted by the square brack- 
ets. To perform the replica average at the level of 
the CW (and not at the level of energy moments) 
is motivated by the general rule that quenched av- 
erages should be performed at the level of the free 
energy and not the partition function. Finally, we 
determined the maximum, C max = C(/3c ma x)> for 
each lattice size and studied the finite-size scal- 
ing (FSS) behaviour of C max and /3c max - The er- 
ror bars on the two quantities entering the FSS 
analysis are estimated by jack-kniving over the 64 
replicas. 

The analysis of the magnetic susceptibility, 
x(/3) = /3N ([(m 2 ) — (m) 2 ]) and the energetic 

Binder parameter V(0) = 1 - [(e 4 )] /3 [(e 2 )] 2 , 
proceeds exactly along the same lines, yielding 
Xmax and /3xmax as wen as ^min and 0v min ■ In or- 
der to be prepared for the possibility of a second- 
order phase transition, the magnetic Binder pa- 
rameter, U((3) = 1 — [(tj 4 )] /3 [(to 2 )] , was also 
measured as, in this case, its crossing for different 
lattice sizes provides an alternative determination 
of the critical coupling j3 c , and the FSS of either 
the maximum slopes or the slopes at (3 C can be 
used to extract the correlation length exponent. 

If the transition was of first-order one would 
expect for large system sizes an asymptotic FSS 
behaviour of the form C max = ac+bcN+. . ., and 
Xmax = o-x + ^x-^ + • • •• However a linear scaling 
with N is not consistent with our data. Further- 
more, at a first-order phase transition one would 
expect that the energetic Binder-parameter min- 
ima approach in the infinite-volume limit a non- 
trivial value related to the latent heat. Our data, 
however, approaches the trivial limit of 2/3, also 
strong evidence for a continuous transition. 
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Figure 1. Maxima of C vs N. The fit line is 
discussed in the text. 



Being thus convinced that the transition is con- 
tinuous, the next goal is to determine the crit- 
ical exponents and the corresponding universal- 
ity class. To this end we have tried to describe 
the scaling of the specific-heat and susceptibility 
maxima with the standard FSS ansatz C max = 
a c + b c N a / Dv and W - a x + b x N^I D \ where 
a, 7, and v are the usual universal critical ex- 
ponents at a continuous phase transition, ac, x 
and bc, x are non-universal amplitudes, and D is 
the intrinsic Hausdorff dimension of the graphs. 
By performing a non-linear three-parameter fit 
to the specific-heat maxima we obtained a/Dv = 
0.22(7), with a reasonable goodness-of-fit param- 
eter Q = 0.10; see Fig. 1. Since the background 
term ac turned out to be consistent with zero, we 
also tried linear two-parameter fits with ac = 
kept fixed and the resulting exponent estimate 
over all data points gave a fully consistent value 
of a/Dv = 0.222(7) (with Q = 0.15) but, as ex- 
pected, a much smaller error bar. The suscepti- 
bility maxima grow very fast with N, such that 
also here the constant term a x can safely be ne- 
glected. The linear fit over all data points yielded 
j/Dv = 0.732(10) (with Q = 0.27), and omitting 
the N = 250 point we obtained 7 / ' Dv = 0.719(14) 
(with Q = 0.31). 

The pseudo-transition points /3c ma x > /^xmax > an d 
/V min were fitted to the standard FSS ansatz 
/3 Cmax = Pc + ccN- 1 ' ", etc. By taking the aver- 
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age of the three estimates of f3 c we estimated the 
transition point to be (3 C — 2.2445(20). This value 
is consistent with the estimate obtained from the 
crossings of the magnetic Binder parameter U 
for different lattice sizes. The estimate of l/Dv 
from the FSS of the pseudo-transition points is 
not very stable but a more precise estimate can 
be obtained by analyzing the FSS of the mag- 
netic Binder-parameter slopes in the vicinity of 
[3 C . Both the maximum slopes and the slopes 
at p c are expected to scale as dU/d/3 cx N 1 / Du . 
The fit to the maximum slopes for N > 500 
gave l/Dv = 0.616(29) (with Q = 0.29), and for 
the slopes at (3 C — 2.2445 we obtained l/Dv = 
0.614(30) (with Q = 0.78). 

On annealed gravity graphs it was found that 
the 10-state Potts model appeared to display the 
exponents of the 4-state Potts model coupled to 
gravity, which are listed as "annealed q = 4" in 
Table 1 below. The values for the exponents mea- 
sured in our simulations appear to be best fitted 
by the "quenched q — 4" exponents derived in 
and also listed in Table 1. On the basis of the 



Table 1 



Critical exponents for Potts models 



Type 


a/Dv 


j/Dv 


l/Dv 


Annealed q = 4 





1/2 


1/2 


Quenched q = 4 


0.177... 


0.709... 


0.589... 


Measured q = 10 


0.22(1) 


0.719(14) 


0.61(3) 



numerical evidence in this paper and the earlier 
results of (|] it would thus seem that the con- 
nectivity disorder of both quenched and annealed 
$ 3 gravity graphs has the effect of softening the 
first-order transition of the q — 10 Potts model 
to a continuous transition. Unlike the case of 
bond disorder, however, we find exponents asso- 
ciated with the 4-state Potts model rather than 
the Ising model. It is perhaps worth remarking 
that in || the 200-state Potts model on annealed 
gravity graphs did apparently display Ising-like 
exponents, so a similar investigation on quenched 
graphs might be interesting. 



3. CONCLUSIONS 

To summarize, we have obtained strong numer- 
ical evidence that due to connectivity disorder 
the phase transition in the 10-state Potts model 
on quenched random gravity graphs is softened 
to a continuous phase transition. This result is 
in contrast to a recent simulation of the 8-state 
model on Poissonian Delaunay/Voronoi random 
lattices where the transition stays first order as on 
regular lattices 0]. It is, however, in qualitative 
agreement with the quenched random bond case 
and like the simulations in || of the 10-state 
Potts model on annealed 2D gravity graphs we 
see exponents associated with the 4-state Potts 
model, in this case the "quenched" exponents. 
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